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Abstract. 

We investigate the dual of the k = Gonihedric Ising model on a 3-D cubic lattice, 
which may be written as an anisotropically coupled Ashkin- Teller model. The original 
K = Gonihedric model has a purely plaqucttc interaction, displays a first order 
transition and possesses a highly degenerate ground state. 

We find that the dual model admits a similar large ground state degeneracy as a 
result of the anisotropic couplings and investigate the coupled mean field equations 
for the model on a single cube. We also carry out Monte Carlo simulations which 
confirm a first order phase transition in the model and suggest that the ground state 
degeneracy persists throughout the low temperature phase. Some exploratory cooling 
simulations also hint at non-trivial dynamical behaviour. 
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1. Introduction 

The Gonihedric Ising model has an interesting history, having originally been formulated 
as a fixed lattice version of a discretized string/triangulated random surface action 
suggested by Savvidy et.al, whose action was given by [l] 

(ij) 

where 0{aij) = \TT—aij \ and aij is the dihedral angle between the embedded neighbouring 
triangles with a common link (ij). The |Xj — Xj\ are the lengths of the embedded 
triangle edges as shown in Fig. ([T]). The aim of this action was to weight the edges of 
non-coplanar adjoining triangles on a discretized surface, rather than the triangle areas 
as is the case with a Gaussian action, in an attempt to search for a continuum limit 
which might be related to a string theory. The discretized random surfaces formed from 
gluing together triangles were intended to model Euclidean string worldsheets. 



X 

1 




Figure 1. Two adjacent triangles in a triangulation of a surface with a common edge 
(ij), showing the co-ordinates of the endpoints Xi j and the dihedral angle ay- 
Translating this action onto a fixed cubic lattice and asking that the requisite 
surfaces be represented by the plaquettes of spin cluster boundaries in some Ising-like 
model trivializes the edge length \Xi — Xj\ dependence, so the statistical weights of 
surface configurations depend solely on the 6{aij) = \t[ — aij\ factors, where the aij are 
now restricted to multiples of 7r/2 radians. The statistical weights of such plaquette 
surface configurations will be determined entirely by the number of bends and self- 
intersections they contain. 

The mapping between the energies of a gas of plaquette surfaces and generalised 
Ising models was studied in some detail by Cappi et.al. [2] who calculated the energies 
of spin cluster boundaries on the 3D cubic lattice for an Ising spin (±1) Hamiltonian 
which contains nearest neighbour (i, j), next to nearest neighbour and plaquette 

k, I] terms 

- (3H = ffiffj + ^2 (yi<yj + O-jO-jCrfcCr^. (2) 

The couplings in such models can be related to the couplings for the area energy of 
plaquettes in spin cluster boundaries, Pa, the energy cost of a right-angled bend between 
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two such adjacent plaquettes, Pc, and the energy cost, /3/, for the intersection of four 
plaquettes having a link in common 

13a = 2Ji + 8J2 
Pc = 2^/3 — 2J2 

/3, = - 4 J2 - 4 J3 . (3) 

The original 3D Gonihedric model [s] constitutes a particular one-parameter slice of this 
family of Hamiltonians: 

K s-^ 1 — K 

2 



H 



-2k CTj 



{{ij)) 



(4) 



For this ratio of couplings /3a = 0, which means that the edges and intersections of spin 
cluster boundaries are weighted rather than their area, which is the antithesis of the 
usual 3D Ising model with only nearest neighbour spin interactions where Pi = Pc = 0. 
The energy of the spin cluster boundaries for the Gonihedric model on a cubic lattice is 
simply given hj E = n2 + Ann^ , where n2 is the number of links where two plaquettes 
on a spin cluster boundary meet at a right angle, n4 is the number of links where four 
plaquettes meet at right angles and k is the free parameter. It is worth remarking that 
the language we have employed implicitly assumes that spin cluster boundaries can be 
clearly identified. As we shall see below this may not be such a simple matter for the 
dual Gonihedric model (it is similarly complicated for the original Gonihedric action, at 
least when n = 0). 

When K = the Gonihedric Hamiltonian becomes a purely plaquette term 



H 



(5) 



which is not the 3D gauge Ising model, since the spins live on the vertices rather than the 
edges of the lattice. This plaquette action displays a first order transition surrounded 
by a region of metastability |4j. It also displays interesting dynamical behaviour with a 
dynamical transition at the lower boundary of the metastable region which appears to 
have many glassy characteristics (5]-[8]. This is intriguing because there is no quenched 
disorder in the Hamiltonian. 

The dual to the k = Gonihedric Ising model was constructed by Savvidy and 
Wegner [9] . They considered a high temperature expansion of the partition function for 
the Hamiltonian in equ.([5]) 

Z(/3)= ^exp(-/3if) 

= 5^ n 

which can be written as 



1 + tanh 



/3 



[OiOjakCri) 



(6) 



Z{P) 



2 cosh 



/3 



3L3 



E 

{s} 



tanh 



/3 



n(5) 



(7) 
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on an cubic lattice. The sum runs over closed surfaces with an even number of 
plaquettes at any vertex and ^(5*) is the number of plaquettes in a given surface. This 
ensemble can be constructed from three differently oriented elementary "matchbox" 
surfaces of the form shown in Fig. (2) along with the unshaded cube. A spin variable 
representing each matchbox then sits at the centre of the cube on the dual lattice and any 
surface in the ensemble can be constructed as a product of the elementary matchboxes 
and unshaded cubes. When two shaded matchbox faces overlap they annihilate to give 
an unshaded face, so the shaded faces can be thought of as carrying a negative sign. 
The low temperature expansion of the dual Hamiltonian 

Hdual = -\Y1 ^'^i ~ ^ 5Z ^^^^ ~ ^ 5Z ^■'■^^ 

(ij) (ik) (jk) 

gives precisely this structure, where a, t and rj represent each of the possible matchbox 
orientations. In equ.(|8]) the spins a, t and r] live on the vertices of the dual lattice and 
the sums are along its orthogonal edges ij,ik and jk. They satisfy 




Figure 2. One of the three possible orientations of an elementary matchbox surface. 



ea = a , er = t , et] = t] 

or^ = = r]^ = e (9) 
ar = rj , tt] = a , rja = t 

with e representing the unshaded cube, which means the spins live in the fourth order 
Abelian group. For convenience in simulations the spins may also be considered as Ising 
(±1) spins if we set rji = ai Ti, which recasts the Hamiltonian into 

Hdual = -^JZ^i^i - ^X^^i^fe - ^^^^j^^kTjTk, (10) 

{ij) {ik) {jk) 



which is recognizable as an anisotropically coupled Ashkin- Teller 10 model with equal 
couplings. Without the four-spin term this would simply be two uncoupled ID Ising 
chains arranged in perpendicular directions and thus display no transition(s), but as we 
see below the coupling via the Ashkin- Teller energy term in the third direction gives 
non-trivial behaviour. 



In the isotropically coupled case the ratio in equ. ( 10 ) corresponds to the increased 
symmetry point of the standard Ashkin- Teller model where the generic Z2 XZ2 symmetry 
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is promoted to Z4. This can be seen explicitly by rewriting the Hamiltonian in terms of 
the four double spins Si = (±1, ±1) to give the 4-state Potts Hamiltonian 

H dual, isotropic = ^ (4^5^,5^ — l) (H) 

where the sum now runs over all the edges orientations. 4-state Potts critical behaviour 
(i.e. a first order phase transition in 3D [TT]) is thus found in the isotropic case. 

In the remainder of the paper we consider the behaviour of the 3D dual Gonihedric 



Hamiltonian, as formulated in equ.(lO) as an anisotropic Ashkin- Teller model, in its 
own right. We first look at the ground state structure of the model, highlighting the 
similarities with the plaquette action and in the light of this discuss coupled mean field 
equations on a cube. We then report on Monte Carlo simulations which are sufficient 
to confirm the nature of the phase transition and look at the effect of different cooling 
rates on the low temperature behaviour. The ground state, mean field and Monte- 
Carlo investigations all highlight the difficult of formulating a standard magnetic order 
parameter, and we discuss the implications. 

2. Ground State 

In the isotropically coupled Ashkin- Teller model with equal positive couplings four 
equivalent magnetized ground states are possible, with the (cr,r) spins taking the 
values (±, ±) at every site, so only paramagnetic or ferromagnetic behaviour is seen 
for the individual cr and r spins in this coupling regime. To investigate the ground 
state/zero-temperature structure of the dual Gonihedric model we use an approach 
which proved useful for the original undualized Gonihedric model and consider possible 
spin configurations on an elementary cube [2] . This is sufficiently large to capture non- 
trivial structure and the full ground state is then obtained by tiling the 3D cubic lattice 
with compatible single cube configurations. 

The full lattice Hamiltonian may be written as a sum over the individual cube 
Hamiltonians he, 

he = -^X^^i^i - - ^X^^i^fc^i^fc, (12) 

(id) {j,k) 

where the additional symmetry factor of | takes account of one edge being shared by 
four cubes. If a configuration of spins on a cube minimizes he the full lattice ground 
state energy density will simply be given by he- 

Looking at the configurations in Fig. ^ immediately makes it clear that there 
is considerably more freedom for possible ground states in the anisotropically coupled 



Hamiltonian of equ.(lO) than in the standard isotropic Ashkin- Teller model. In addition 
to the reference ferromagnetic ground states in Fig. ([3^) it is possible to fiip a single 
horizontal face of r spins on the cube at zero energy cost as in Fig. (|3}d), or a vertical 
face of either a (Fig. ^)) or both a and r spins (Fig. ([3]i)) on the differently oriented 
vertical faces. Flipping two faces in the same orientation takes one between the different 
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possible ground states of the standard isotropic model, so it is the ability to flip a single 
zero energy cost which confers the extra freedom in the anisotropic model. 

It is also possible to combine the differently oriented single face flips on the cube 
without increasing the energy of the conflguration. Tiling the entire lattice with such 
conflgurations then shows that ground states may contain flipped planes of cr, r or 
err spins (depending on the orientation) with respect to reference purely ferromagnetic 
ground states. It is possible for two orientations of flipped spin planes to intersect 
pairwise along a line or for three differently oriented planes to intersect at a point. The 
distribution of flipped spin planes in a ground state is thus completely arbitrary. 

The ground state degeneracy of the dual Gonihedric model is thus similar to 
that of the original plaquette Hamiltonian of equ. ([S]) where (possibly orthogonal, 
intersecting) planes of spins may also be flipped at zero energy cost. Intriguingly, this 



++ 



(a) 



+ 1- 




++ 




ax 



Figure 3. Some possible ground state spin configurations on a cube, the cr, r values are 
shown at each site. The directions of the couplings in the Hamiltonian are indicated, 
as are the faces on which spins are flipped. 



leads to the same difficulties in defining a suitable magnetic order parameter for the 
dual Hamiltonian as one faces with the plaquette Hamiltonian. Since arbitrary, and 
arbitrarily separated, spin planes may be fiipped at zero energy cost a standard, or even 
a staggered magnetization, will generically be zero in such a state. 

Of course, it is not guaranteed that ground state/zero-temperature degeneracies 
are maintained at finite temperature. Indeed, for the original Gonihedric model low 



temperature expansions by Pietig and Wegner 12 showed that when k ^ 0, where 



it was still possible to fiip arbitrary parallel spin planes to give a sandwich ground 
state, the ferromagnetic state had a lower free energy at finite temperature. The 
degeneracy did, however, persist at finite temperatures for the k = model. In the 
case of the dual Gonihedric model the Monte-Carlo simulations discussed below find 
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(cr) = (r) = ((Tt) ~ in the low temperature phase, though there is a clear sign of 
a phase transition in the energetic observablcs and also the magnetic susceptibilities. 
This suggests that the flip symmetry persists at finite temperature throughout the low 
temperature phase. 



3. Mean Field 



In the mean field approximation an expression for the free energy is written by replacing 
the exact values of the spins with average site magnetizations and adding an entropy 
term. To take account of non-trivial structure in such a calculation one can again work 
at the level of the cubes, as for the ground state, and write down the coupled equations 
for the sixteen site magnetizations on the cube. The calculation of the mean field free 
energy in this manner is thus a direct elaboration of the method used to investigate the 
ground states. The total mean field free energy is written as a sum of the elementary 
cube free energies (f){'mc,nc), given by 

/3(j){mc,nc) ^ ~ f XI ^i^J ~ ^ riiUk - ^ mjmkrijnk 

{i,j)<ZC {i,k)cC {j,k)GC 

+ ^ + rrii) ln(l + rm) + (1 - rm) ln(l - m^)] 

+ ^ + rii) ln(l + m) + (1 - m) In (1 - m) ] (13) 

icc 

where mc, ric is the set of magnetizations of the elementary cube with mj, rii the average 
site magnetizations for the Ui and Tj spins respectively. The log terms give the entropy 
for each of the types of spin. Minimizing this free energy gives a set of sixteen coupled 
mean-field equations 

d(f){mc, nc) 



= 



d<Hrnc^ = (14) 

orii (i=i...8) 

(one for each corner of the cube and spin type) rather than the familiar single mean 
field equation for the standard nearest neighbour Ising action. The resulting equations 
are all of the form 

mi = tanh[/3(m4 + m2 ni 712)] 

ms = tanh[^(m5 + myny^s)] (15) 
rii — tanh[^(n5 + 712 mi m2)] 

Us = tanh[/3(n4 + nymymg)] 
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where we have labelled the magnetizations on a face of the cube counter-clockwise 
1 . . . 4 and similarly for the opposing face 5 ... 8, as shown in Fig. Q. If we solve these 
equations iteratively at different temperatures we arrive at zeroes for a paramagnetic 
phase or various combinations of ±1 for the magnetized phases on the eight cube vertices. 



5 , 



Figure 4. The labelling of sites used in writing the mean field equations for the cube. 
The directions of the mean field spin couplings in the Hamiltonian are again indicated. 



A potential problem with an iterative scheme to solve the system of mean field 
equations 

mf+^) = f,[m', n'^] , = f,[m', n'^] , (16) 

is that it might fail to converge if an eigenvalue of drnf^^^^ /dm^^ or dnf^'^^^ /drij is less 
than — 1 j2]. This is easily remedied by modifying the equations to 

(fc+i) _ {Mm\n''] +amf) 



n. 



l + a 
l + a 



(17) 



for a suitable a, and we have employed this here to ensure stability. 

The mean field solution is then given by gluing together the elementary cubes 
consistently to tile the complete lattice, in the manner of the ground state discussion. 



In the limit /3 — > oo the mean field equations of equ.(16) become the system of equations 
nil = sgn[m4 + m2 rii 



ms = sgn[m5 + rrij nj ns] 



(18) 



rii = sgn[n5 + n2 mi 7712] 



77-8 = Sgn[n4 + 7717 ^s] 
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which are, as they should be, compatible with the various ground state structures shown 
in Fig. ([3]). Solving the mean field equations numerically finds a (single) transition at 
13 ~ 0.83 from a paramagnetic state to one of the possible ground states. If the iteration 
is seeded with spin values close to ±1 one of the ferromagnetic states is picked. 

It would be interesting to refine the mean field solution further by employing 
a cluster variational approximation, which essentially amounts to "improving" the 
entropy term and can be combined with Fade approximant methods to obtain quite 
accurate critical exponent estimates 13 . This has been done successfully for the original 
Gonihedric model 14 , but we do not pursue this further here, turning instead to Monte 
Carlo simulations to sketch out the phase diagram of the model. 



4. Some (Modest) Monte Carlo 



For comparison purposes the phase diagram for positive couplings for an isotropically 
coupled Ashkin- Teller model is shown schematically in Fig.(|5]) where J2 is the two- 
spin and J4 is the four-spin coupling 10,11 . Decreasing the temperature along the 



J2 = Ji line moves from a paramagnetic phase to a phase in which all of (a), (r) 
and (ar) ("polarization") are non-zero, which is sometimes called the Baxter phase. 
The transition takes place at the four state Fotts point and is thus first order, as 



we saw from a direct rewriting of the Hamiltonian in equ.(ll). We are principally 



0.4 - 



<a> <T> <ax> 



<ax> 



0.4 



Figure 5. A schematic drawing of tlie pliase diagram of the isotropic 3D Ashkin- Teller 
model for positive two-spin, J2, and four-spin, J4, couplings. The 4-state Potts point 
on the J2 = J4 line is marked as P. The indicated order parameters are non-zero in 
the phases shown on the diagram, and first order transition lines are shown as dashed 
lines, second order lines as solid. The paramagnetic phase surrounds the origin and the 
effect of decreasing the temperature along the J2 = J4 line is indicated by an arrow. 



interested in determining the nature of the transition in the dual Gonihedric model 
here, for comparison with the phase diagram of the isotropic Ashkin- Teller model in 
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Fig. ^ rather than carrying out a high accuracy scaling analysis, so we use relatively 
modest lattice sizes and statistics in our simulations and employ a simple Metropolis 
update. Lattices of size L = 10, 12, 14, 16, 18 and 20 with periodic boundary conditions 
for both the a and r spins were simulated using both hot and cold starts at various 
temperatures. Following a suitable number of thermalization sweeps determined by the 
energy autocorrelation time, 10^ measurement sweeps were carried out at each lattice 
size for each temperature simulated. 

An estimate for the phase transition point of the original plaquette k = 



Gonihedric model is the value in 15 which takes account of the (effectively) fixed 
boundary conditions employed in the simulations there to fit to a suitable scaling 
form with the correct leading and subleading finite size corrections in such a case, 
f3^{L) = I3c + ai/L + as/L^. This found 13^ = 0.54757(63). Allowing for factors of 2 
in the coupling definitions, an estimate for the dual transition temperature /3* is then 
given by the standard formula (3* = — ln[tanh(/3c/2)] = 1.32. 

A plot of the energy is shown for various lattice sizes in Fig. ^ where there is 
clearly a sharp drop in the region of /3 ~ 1.38, somewhat higher than the estimate from 
the dual transition temperature. We can get a rough idea of the energy autocorrelation 




1 .22 1 .24 



Figure 6. The energy for variously sized lattices ranging from 10"^ to 20'^ from left to 
right. The lines joining the data points are drawn to guide the eye and hot starts have 
been used in all the simulations. 



time Te in the simulations by comparing the naive estimate for the variance (i.e. specific 
heat) 



{E, - {E)y 



(19) 



1 

where is the number of measurements carried out, with a jack-knifed estimate using 
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binned data ejK- The two are related by 

^JK = \l -(-naive- (20) 

V nm 

Away from the transition point we find Tg ~ 1 but large values of ~ 10^ appear in its 
vicinity. 

There is no signal of the observed phase transition in any of the standard magnetic 
order parameters (a), (r) and (cxr). This suggests that the degeneracy observed in 
the ground state features at finite temperatures as well. The susceptibility for both 
the individual a and r spins and the polarization ar is, however, non-zero and does 
show a signal at the phase transition point where, like the energy, it drops sharply. In 
Fig. ([7]) we plot the polarization susceptibility for 10^ — 20^ lattices. The behaviour of 
the individual spin susceptibilities is very similar, with that for both the a and r spins 
showing a sharp drop from the high temperature phase to a much lower value at the 
same point. 











10 \ 12 ' 


Si 4 
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18 


20 













1.26 1.28 1.3 1.32 1.34 1.36 1.38 1.4 



Figure 7. The polarization susceptibility, x for a 10"^ to 20"^ lattices, showing a sharp 
drop near the pseudo-critical point in the region oi (i — 1.375 for the various lattice 
sizes. As for the energy in Fig. ^ the lines between the data points are drawn to 
guide the eye and hot starts have been used. 



A good indicator of a first order transition is a bi-modal energy distribution at the 
transition point, so we also histogrammed the energy during the simulations. Looking 
at an energy histogram from a simulation sufficiently near the finite size pseudo-critical 
temperature should display a two-peak structure for a first order transition. A typical 
example for a 10^ lattice close to its finite size pseudo-critical temperature at [3c = 1.275 
is shown in Fig. ([s]), where we have histogrammed the 10'' measurements of the energy 
which were carried out after each full lattice sweep of the a and r spins. The double peak 
structure expected of a first order transition is clearly visible. A direct consequence of 
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Figure 8. The energy histogram from a simulation with 10'' sweeps on a 10"^ lattice 
near the finite size transition point (in this case /3 = 1.275). P{E) is shown on a 
logarithmic scale. 



the bi-modal energy distribution near criticality is a non-trivial limit for Binder's energy 
cumulant. This is defined as 

which approaches 2/3 at a second order transition point and a non-trivial limit at a first 
order point. We observe a non-trivial minimum value of Ue of 0.60(2), which varies 
little across the lattice sizes simulated. 

Similarly, the (3 value of the minimum of C/e on an lattice, (3rnin{L), is expected 
to scale as /3min{L) = Pc — 0{l/L^) for a first order transition. If we plot the estimated 
minima positions for the various lattice sizes against 1/L^ we get a reasonable fit to this 
behaviour with a value of Pc ~ 1.388(4) and a Xdof of 1-34 when the smallest lattice 
size of 10^ is dropped from the fits. The estimated value of /3c from this procedure is 




Figure 9. The scaling of the position of the minimum of the Binder energy cumulant 
against the inverse volume. 
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consistent with the behaviour of the energy and susceptibility jumps but, as we have 
already noted, it is somewhat higher than the dualized value calculated from from the 
measurements in 



15 



There is strong hysteresis around the transition point which may affect the accuracy 
of any such estimates. We show the result of using both ordered (cold) and disordered 



(hot) starting configurations for a 14^ lattice in Fig. (10) with, in both cases, relaxation 



times of 10^ sweeps followed by 10^ measurement sweeps. This behaviour is again 











COLD 


HOT 









1.24 1.26 1.28 1.3 1.32 1.34 



Figure 10. Measurements of the energy for hot and cold starts in the region of the 
transition point on a 14^ lattice. 

strikingly similar to that seen in the original plaquette Hamiltonian js]. 
5. Some (Very Modest) Dynamics 

An intriguing feature of the original n = Gonihedric model is its highly non-trivial 
non-equilibrium behaviour, including a region of metastability around its first order 
transition point and a dynamical transition which displays many glassy characteristics. 
We have seen that the dual Gonihedric model's equilibrium phase diagram is similar to 
the original, so it is natural to inquire whether this similarity also holds for dynamical 
behaviour. 

As a start in this direction, we cooled differently sized lattices at various rates 
from disordered (hot, "liquid") starts in order to see if there was any evidence of the 
potentially glassy behaviour seen with the plaquette Hamiltonian. We cooled lattices of 
size 20^, 40^ and 60^ starting with a disordered configuration at a temperature of T = 3. 
From Figs. ( |Tl|l2| it is clear that there is little difference between the 40^ and 60^ lattice 
results, while the 20'^ lattices may still be subject to stronger finite size effects. With 
regards to the numerical estimates of critical temperatures, this also suggests that the 
equilibrium simulations in the previous section may have been carried out on rather 
small lattices. 
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In Fig. (11) we can see that with a slow cooling rate of 5T = 0.00001 per sweep, the 



systems still relax to a (ground) state with E = —1.5. The jump in the energy at the 
phase transition seen in the time series at T ~ 0.72 on the larger lattice sizes is consistent 
with the estimate of /3c = 1.388(4) obtained from extrapolating the Binder cumulant 
values. Fig. (12), where a faster coohng rate of 5T = 0.001 per sweep is employed, is 




Figure 11. The time series of energy measurements obtained from cooling 20'^, 40^ 
and 60'^ lattices from a hot start at a rate of 6T = 0.00001 per sweep. 



perhaps more interesting. Once again there is little difference between the 40^ and 60^ 
lattices, but this time they do not relax to the ground state energy of £^ = —1.5, but 
are trapped at a higher value. Similar results in the plaquette model IsIIg] were taken 



Figure 12. The time series of energy measurements obtained from cooling 20'^, 40^ 
and 60'^ lattices from a hot start at a rate of 6T = 0.001 per sweep. 



as an indicator of possible glassy behaviour, though there is evidence that what is seen 
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is an echo of a mean-field spinodal point |8] at which the supercooled high-temperature 
("liquid") phase becomes physically irrelevant. 

6. Discussion 

We have studied the dual of the k = Gonihedric Ising model in three dimensions, 
which may be formulated as an anisotropically coupled Ashkin- Teller model. We noted 
that this dual Gonihedric model displays clear signals of a first order transition, such 
as a bimodal energy histogram and a non-trivial limit for Binder's energy cumulant 
(like the isotropically coupled equivalent), but it has a highly degenerate ground state 
(txn-like the isotropically coupled equivalent). 

The magnetic order parameters (a), (r), (ar) show no strong signal at the 
phase transition point which is, however, clearly visible in the energy and various 
susceptibilities as well as the Binder (energy) cumulant. The absence of conventional 
ferromagnetic order suggests that the degenerate ground state structure persists to finite 
temperatures, as is the case for the k = Gonihedric model. It would be an interesting 
exercise to carry out a low temperature expansion of the dual Hamiltonian to verify 
this by comparing the energy of the states with fiipped spin planes to a ferromagnetic 



reference state in the manner of 12 



Another aspect which merits investigation is the crossover to the isotropic model, 
which has a much simpler ground state structure and known, simple order parameters 
for the low temperature phase(s). In any such endeavours the form of the anisotropic 
action suggests that it might be more amenable to a cluster simulation than the original 
plaquette action. The first order nature of the transition means that this would not 
confer such great advantages over local updates as in the case of continuous transitions, 
though cluster updates could be employed in conjunction with multihistogramming 
methods of various sorts for maximum numerical efficiency. 

To investigate the non-equilibrium behaviour of the model, however. Metropolis (or 
other local) dynamics should be employed. We have made a start in this by conducting 
some cooling experiments which show that the phenomenology of the dual Hamiltonian 
appears to be remarkably similar to that of the original k = plaquette Hamiltonian. 
Under very slow cooling the ground state energy is achieved, but faster cooling appears 
to trap the system in a higher energy state. More extensive simulations along the lines of 
those conducted in [7||8] for the plaquette model and the the coupled two layer system 
(CTLS) would be useful to discern whether the "bubbling and coarsening" scenario 
posited there for the low temperature behaviour also applied in the case of the dual 
Gonihedric model, or whether more conventional coarsening dynamics was seen. 

If the low temperature behaviour does, indeed, display (pseudo-)glassy 
characteristics it would also be useful to elucidate the nature of the self-induced 
frustration which is presumably causing it. For both the plaquette Gonihedric model and 
the CTLS multi-spin interactions appear to play a vital role. In view of the simplicity 
of the Hamiltonian, the dual Gonihedric model might also provide a further test case 
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in which to explore the approach of 16 , which links dynamical, glassy behaviour in 



classical systems to quantum phase transitions. 
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